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^ (54) Title: METHOD AND SYSTEM FOR BALANCING THRUST DEMANDS 

(57) Abstract: In the control of spacecraft, thrusters must generate a commanded force or torque that is computed from the control 
algorithms in order to maintain a certain attitude or orbit In general, there is an infinite number of solutions as to how the thrusters 
are fired in order to achieve the commanded value. However, additional imposed constraints on the mass flow allow only a small 
subset of those solutions to be feasible. By way of example, two constraints are considered in this invention: first, the mass flow of 

— each thruster must exceed a minimum value and second, the total mass flow of all thrusters for any generated force/torque must be 

J2 constant If only a non-feasible solution is found, known methods turn it into a feasible one by reducing the commanded force/torque 
until the constraints are met Reducing the commanded force/torque is, however, a disadvantage because the performance may 

^5 decrease and even lead to instability of the closed loop system. This invention concerns a simple and time efficient method/system 
for improving an unfeasible solution in such a way that reducing the commanded force/torque can be minimised or even completely 

^ be avoided and, therefore, higher feasible generated forces/torques than those with known methods can be generated. 
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METHOD AND SYSTEM FOR BALANCING THRUST DEMANDS 

FIELD OF THE INVENTION 

The present invention relates to a method and system for balancing 
5 thrust demands, and more particularly, but not exclusively, concerns a method 
and system for balancing thrust demands for attitude and translation control of 
spacecraft. 

BACKGROUND OF THE INVENTION 

An important issue during the design of an attitude and translational 
10 control system for spacecraft is the modulation of commanded control forces 
and torques towards the thrusters. The design gets even more challenging if 
additional constraints are taken into account, for instance the requirement that 
the total mass flow remains constant and that each thmster generates a 
predefined minimum thrust. 

15 The known modulating algorithm proposed in Hai Ping Jin. Peter Wiktor, 

and Daniel B Debra's report, "An Optimal Thruster Configuration Design and 
Evaluation for Quick Step," IS"* IFAC Symposium - Automatic Control In 
Aerospace - Aerospace Control 1994, September 12-16 1994 aims to tackle 
the design issue mentioned above. A principal drawback of this known method, 

20 however, is that the constraints cannot be satisfied for high control demands 
with respect to the ideal maximal thrust which can be generated from the 
thrusters. The reason for this drawback is that only a small fraction of the 
thmster null-space is used for the modulation. In contrast, the null space 
search algorithm of the present Invention, as proposed hereinafter, uses the 

25 complete null space and therefore higher control demands while satisfying 
additional constraints can be generated. Another advantage associated with 
the present invention is that the proposed null space search algorithm is fast 
and compatible for application to any modulator and therefore can be used to 
improve any thruster signal proposed by any modulator. 

30 
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nR.IFnXR AND SUMMARY OF TH E INVENTION 

The present invention aims to overcome or at least substantially reduce 
some of the above-mentioned drawbacks. 

It is an object of the present invention to provide an improved method 
5 and system for balancing thmst demands with a capability to enhance the 
performance for thruster modulation by generating higher feasible 
forces/torques than would othenwise be generated using known 
methods/systems, whilst satisfying the prescribed constraints. 

It is another object of the present invention to provide a simple, cost- 
10 effective and time efficient method of utilising the complete null space, enabling 
high control demands to be generated whilst satisfying additional constraints. In 
this connection, it is to be understood that the null space search algorithm for 
use in the present invention is fast, reliable and compatible for application to 
any modulator and can be effectively used to Improve upon any thruster signal 
15 proposed by any modulator. 

In broad terms, the present invention resides in the concept of generating 
a number of thrust command signals for a plurality of thmsters and processing 
the thrust command signals in a simple, logical fashion so as to (1) identify 
exacUy which of the thmst demands associated with the thmsters do or do not 
20 satisfy predetermined constraints and (2) balance the various thrust demands in 
such a way that each of the thrusters can operate in accordance with the 
predetermined constraints. 

Therefore, according to the present invention, there is provided a method 
of balancing thrust demands comprising the steps of: (a) controllably generating 
25 a number of thrust command signals for application to a plurality of thmsters; (b) 
comparing the levels of the generated thmst command signals in relation to a 
number of predetemiined constraints; (c) identifying which of the thmst 
demands associated with said plurality of thmsters can operate in accordance 
with said predetermined constraints and which of the thmst demands 
30 associated with said plurality of thmsters can not operate in accordance with 
said predetemiined constraints; and (d) balancing the various thmst demands 
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such as to permit each of the thmsters to operate in accordance with said 
predetermined constraints. 

In this way, available thrust is used in the invention to a higher extent 
than in conventional techniques and at the same time, the prescribed 
5 constraints of the thrusters are suitably satisfied. 

In accordance \Arfth an exemplary embodiment of the invention which will 
be described hereinafter in detail, the balancing step is perfonned by shifting a 
predetermined amount of thrust demand associated with particular thrusters 
operating in accordance with said predetermined constraints to one or more of 
10 the thrusters which do not operate in accordance with said predetermined 
constraints. 

Advantageously, the balancing step is effected by application of a 
predetermined combination of null space vectors. 

Advantageously, the predetermined constraints are selected so that (i) 
15 the thrust demand associated with each of said plurality of thrusters is more 
than a predetermined value and (ii) the total mass flow associated with said 
plurality of thrusters is constant. 

Conveniently, the method of the invention further comprises the step of 
modulating a number of forces/torques to provide said number of thrust 
20 command signals. Further, it is to be appreciated that the balancing step and 
the modulation step are conveniently effected independently of each other. 

Advantageously, the modulation step is effected by means of a non- 
linear pseudo-inverse modulator. 

The present invention extends to a system for balancing thrust demands 
25 comprising: means for controllably generating a number of thrust command 
signals for application to a plurality of thrusters; means for comparing the levels 
of the generated thrust command signals in relation to a number of 
predetermined constraints; means for identifying which of the thrust demands 
associated with said plurality of thrusters can operate in accordance with said 
30 predetermined constraints and which of the thrust demands associated with 
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said plurality of thrusters can not operate in accordance with said 
predetermined constraints; and means for balancing the various thrust demands 
such as to permit each of the thrusters to operate in accordance with said 
predetermined constraints. 

5 The present Invention further extends to a spacecraft system adapted 

and arranged to carry out the above described method. 

The present invention is advantageously embodied in software. 

The above and further features of the Invention are set forth with 
particularity in the appended claims and will be described hereinafter with 
1 0 reference to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS. 

Figure 1 shows a graphical display of null space vector components for 
use in the Invention; 

Figures 2 to 5 show typical force/torque demand profiles which can be 
15 realised by use of the invention; 

Figure 6 shows a typical error signal force demand - thrust profile which 
can be realised by use of the invention; 

Figure 7 shows a typical thruster time history profile which can be 
realised by use of the Invention; and 

20 Figures 8 to 11 show some typical results for various kinds of modulator 

as obtained with/without the inventive null space search of the invention. 

DETAILED DESCRIPTION OF AN EXEMPLARY EMBODIMENT 

1. THRUSTER MODULATION ALGORITHM 

Having regard to the foregoing, the present invention deals with the 
25 following technical problem: For example, given 6 values of the forceAorque 
demand of the controller, how do the associated 16 thrusters have to be fired to 
satisfy this demand, while at the same time ensuring that the total mass flow 
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remains constant throughout the mission and that each thruster undergoes a 
certain minimum of mass flow? 

The known modulator is the known non-linear pseudo-inverse algorithm 
given in Hai Ping Jin, Peter Wiktor, and Daniel B Debra's report, "An Optimal 
5 Thruster Configuration Design and Evaluation for Quick Step." 13*^ IFAC 
Symposium - Automatic Control in Aerospace - Aerospace Control 1994. 
September 12-16 1994. Further analysis which will be presented in the 
description hereinafter has shown that the performance of this known modulator 
is good unless the force and torque demands are too "high", but it is not able to 
10 use the total thrust available generated by the helium boiling off. It has been 
found that the known modulator fails already at a low force/torque demand level 
whereas the physical thrust available predicts that three times the level could be 
used. Since using the total thrust available would mean that the use of magnetic 
torque rods (MTQ) is no longer mandatory, emphasis is placed hereinafter on 
15 the design of an appropriate modulator. 

In the following description, section 1.1 poses the technical problem to be 
addressed by the present invention. In section 1.2, different kinds of modulator 
are presented for use in. the invention to solve the technical problem. All 
modulators presented hereinafter do not consider directly the prescribed mass 
20 flow constraints. Therefore, there might be zero or even negative mass flow 
components. This issue is taken care of equally for all modulators in sections 
1.3 and 1.4 hereinafter. In section 1.5 it is to be noted that typical force and 
torque profiles are presented which are used for the modulators as open loop 
input signal. Finally, section 1.6 concludes with an evaluation of the generated 
25 results. 

1.1 TECHNICAL PROBLEM STATEMENT 

Find a modulator algorithm that solves the equation 



'COG 



:=b 



6x1 
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(1.1.1) 



'16x1 
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where 

f,t are force/torque demand, 

mfh vi. n are mass flow, direction and location of thruster /, 
g is the gravitational acceleration on Earth surface (-9.81 m/s^), 
ISP=130 s is the specific impulse of the thrusters and 

rCOG is the location of the spacecraft COG w.r.t. the centre of the thruster 
configuration. 

and additionally solves the mass flow equation 

'f,rnf, = 2:c, = wi/ = 1.5 mg/s C -I -2) 

as well as the minimum mass flow constraints 

X 0.009375 mgls<= mf, for i = 1,...16 C -I -3) 

Note that in this example, there are 16 thrusters (hence. 1=1, ...16) but the 
number of thrusters in the invention is not critical so lonq as there are at least 
two thrusters. 
1.2 MODULATORS 

1.2.1 LINEAR OPTIMISATION M ODULATOR 

Several modulators have been developed in this study. In order to 
evaluate the perfomiance of the modulators with respect to the computational 
effort they are compared in their perfonnance with the result of the solution of 
the Linear Optimisation problem; here. Matlab function "linprog" was 
conveniently used to solve the problem 

a^^izeZx, while solving = for Xj^^x, ^x, (1-2.1) 

The description of Matlab function "linprog" claims to find the best 
possible solution for this problem. Therefore, the perfomiance of this modulator 
is used as reference for the other modulators telling us what is physically 
possible to achieve. 
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The reason why the Linear Optimisation modulator was not implemented 
in the STEP-simulation model was the very large computation time, and 
therefore it seems unlikely to use that during the mission. 
1.2.2 NQN-LINEAR PSEUDO-INVERSE MO DUI-ATOR 

5 This thruster modulation algorithm is the known non-linear pseudo- 

Inverse algorithm given in Hai Ping Jin, Peter Wiktor, and Daniel B Debra's 
report, '^An Optimal Thruster Configuration Design and Evaluation for Quick 
Step," 13**' IFAC Symposium - Automatic Control in Aerospace - Aerospace 
Control 1994, September 12-16 1994. This algorithm modulates the thrust 

10 levels of the thrusters to provide desired force and torque vectors, whilst giving 
a constant mass flow rate. 

It has already been found easy In simulations to saturate the thrusters, 
given their very low thrust levels. Typically the 10% of mean flow rate constraint 
is violated (rather than the high thrust limit). As originally implemented, the 

15 modulation algorithm 'cheats' under these circumstances, by providing the 
required force and torque vector, but doing so by exceeding the total flow rate. 
Under very noisy conditions, the flow rate can be much higher than can be 
supplied by normal He boil off. 

1.3 SATISFACTION OF MASS Fi OW EQUAL ITY CONSTRAINT 

20 The mass flow equality (1.1.2) can be enforced, simply by adding 

components of the null space vector: The total mass mtot is scaled by adding to 
the mass flow vector x the vector Xdm^o yield 

where 

25 1 1], (1-3.2) 

then the total mass flow m' is 

;„•= V (x(0 + x,,jr (0) = f.x(i) +f^x^ (0 = m,, + (m/ -m^,) = mf (1 .3.3) 
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aiways satisfying tlie mass flow equality. 

1.4 SATISFACTION OF MINIMUM MASS FLOW CONSTIRAINTS 

The problem for all modulator algorithms is the satisfaction of the 
minimum mass flow constraints. Experiments show that those get easily 
5 violated if the force/torque demand increases. 

1.4.1 EXPLOITING THE COMPLETE NULL SPACE OF THRUSTER 
CONFIGURATION 

So far, the complete null space has not been exploited by any known 
method presented so far. The idea for use in the present invention is now to use 
10 "suitable" combinations of all null space vectors in order to put more thrust on 
those thrusters which violate the prescribed constraint(s) without changing the 
generated forces and torques. 

Mathematically speaking, there is a vector Xf with 

Z:x.(0 = -»/ (1-4.1) 
15 solving 

b = Ax, (1.4.2) 
having at least one component y that violates the minimum mass flow constraint. 

The problem is how to find a specific null space vector Axi with 2 /ixi(})=0 such 
20 that 

x,0)<^0)+Ax,0) (1-4.4) 
and "hopefully" also 

Jc,a) + Ax,0)>x (1.4.5) 
because then the minimum mass flow constraint is satisfied. 
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Say the Matrix A in (1.4.2) has 6 rows and 16 columns and has full rank. 
Then the null space matrix Xo has the dimension 16x10, i.e. any linear column 
vector combination of Xo gives 

A [A, Xo, + Jljo JCo,o ] = 0 (1 .4.6) 

where X■l,...X^Q are any scalars. 

The solution procedure requires three steps. For simplicity but without 
limitation of the method, only one violation is considered. 

Step 1.) Identify the violating thmster, say thruster j. 

Step 2.) In the best case there is a null vector which strongly points in the 
direction of the thruster which violates the constraint. Adding this vector to the 
thmst would then ease the constraint violation. A good guess is the null vector 
Xo which is computed below: 



J 
..7+1 

16 



(1.4.7) 



xl=XoX (1.4.8) 

Note that xl'is a vector pointing "as strongly as possible" towards, 
because A* is the least squares en-or solution of the "ideal" null space vector in 
(1.4.7). Now, from Xo is a vector Ax, derived being properly scaled such that the 
sum of its components is zero: 

Ax,(/c) = k: Ax" (1.4.9) 

T 



J 16 



(1.4.10) 
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Note that k is a free parameter which will be chosen properly. 

Step 3.) Reorder now vector Xi such that the elements decrease from the top to 
the bottom (so the index j becomes now index 16) components and apply the 
same ordering rule to Ax^. The result is illustrated as a graphical display of 
5 vector components in Figure 1 . 

As shown in Figure 1 , it is to be noted that the components of vector xi 
can now be interpreted as discrete function values. The "function" Xi cannot 

increase and the Integral of the "function" Ax^ is always zero. "Hopefully", there 
will be a bulge Ax^ at index 16, such that the violating thrust may get more 
10 thrust from the null space vector than the other components of x^ The technical 
problem addressed by the present invention can now be expressed in a more 
transparent way, namely: 

Find that k that maximises the minimum of the function Xi + k Ax^ . 

To achieve that, the derivative of Xi + kAj^ is computed and set to zero, thus 
1 5 yielding 1 5 possible values for k: 

l-J = 0=>1C = =ir- (1.4.11) 

d'' indices of vectors^' dAx, 

All possible values for k in (1.4.1 1) are computed and this k yielding the largest 
minimum of x^ + KAxj^ is taken. The process is repeated, until there is no further 
improvement. 

20 It turns out that the null space search algorithm for use in the invention 

converges very fast, either to a "positive" solution or that it finds that there is no 
possible improvement. 

The null space search algorithm can be advantageously used on any modulator 
for use in the invention in order to see if constraint violations can be removed. 

25 1.4.2 REDUCTION OF FORCE/TORQUE DEMAND IN CASE OF 
CONSTRAINT VIOLATIONS 
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If the minimum thaist constrained is still violated, the force/torque 
demand will be scaled down until this is no longer the case.. Care must be taken 
such that the reduction Is not larger than necessary. 

In order that the system should behave realistically under such 
5 circumstances, a low thrust saturation diaracteristic has been designed and 
implemented for use in the Invention. The null space search algorithm 
conveniently preserves the direction of both the force and torque demand 
vectors, whilst scaling down their magnitudes to meet the mass flow constraint. 
A similar satufBtion scaling has been used on previous satellite programmes; 
10 this design means that a large demand on one axis does not cause the 
actuation on the other axes to be of incorrect sign, which advantageously 
improves overall robustness. 

1-4.3 DESCRIPTION OF THE INVENTIVE ALGORITHM 

Say the thrust vector xi has a smallest component smaller than the 
15 minimum thrust. The following algorithm for use in the invention scales down 
the demand until the smallest thrust component is hitting just the minimum 
constraint: 

1 .Find index j such with m\r){xi) = XiQ) < x 
2. compute X2 = aj xi + Pj Xo where 

20 "^-^^ (1.4.12) 

where xo is a null space vector. In the STEP study, the following vector was 
used: 

^ : (1.4.13) 
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3.Set now Xi := X2 and goto step 1, else goto 4. 

4.Finished. The scaling down factor of the thmst demand is 

Xd = a<^'a<2\...a^'"> (1.4.14) 

where m is the number of loops from step 1 . until step 3; it can be shown that m 
5 is always less than or equal to 16. 

Having regard to the foregoing, thus, in operation of the method of the 
present invention, the following steps are required to be carried out: 

(1) a number of thrust command signals are controllably generated for 
application to the thrusters (two or more thrusters), 

10 (2) comparing the levels of the generated thrust command signals in relation 
to a number of predetermined constraints, the constraints being of a kind 
as set out hereinabove for example, 

(3) identifying which of the thrust demands associated with the thrusters can 
operate in accordance with the constraints and which of the thrust 

15 demands associated with the thrusters can not operate in accordance 

with the constraints, and 

(4) balancing the various thrust demands such as to permit each of the 
thrusters to operate in accordance with the constraints. 

In this connection, it is to be understood that the balancing step could be 
20 effected as part of the thrust modulation step or alternatively, if desired, it 

could be effected independently of the thrust modulation step. 

1.5 FORCE/TORQUE PROFILE 

In order to test and compare the applicability of the abovementioned 
modulators (see section 1.2) for use in the Invention, four different 
25 force/torque "open loop" demand profiles are considered by way of 

example. 
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1.5.1 PROFILE #1 

Referring to Figure 2, there is sliown a typical STEP demand profile for a 
differential acceferometer. It is extracted fronn closed loop simulation 
runs with all available disturbances and sensor noises switched on. The 
5 torques do not include magnetic disturbances from spacecraft dipole. 

1.5,2. PROFILE #2 

Refem'ng next to Figure 3, there is shown another profile which is the 
same as profile #1 but additionally it does include magnetic disturbances 
from spacecraft dipole. 

10 For the dipole, in each axis +/-5Am^ is assumed. The first parts of the 

profile #2-torques are taken from profile #1, and additionally with the 
Earth magnetic field vector in body co-ordinates extracted from the same 
simulation njn as in profile #1. 8 different magnetic torque sets 
corresponding to all 8 possible dipole component sign combinations are 

15 generated which are subtracted from the torques of profile #1. This 

results in eight sets of forces and torques, where the forces of each set 
are the same and the torques differ in magnetic disturbance torque. 

1.5.3 PROFILE #3 

Referring next to Figure 4, there is shown another typical force/torque 
20 demand which can be realised in the invention. The purpose of this 

profile is to test the modulators under "extreme" conditions. This profile 
is generated simply by setting one thruster each to the maximal 
allowable thrust and the remaining thrust is equally distributed among the 
remaining thrusters. This results in 16 different single force/torque 
25 demands. 

1.5.4 PROFILE #4 

Referring next to Figure 5, there is shown another typical force/torque 
demand which can be realised in the invention. This profile is a very 
special case where all components of the force and torque vector are 
30 nonnegative. 
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1.6 RESULTS AND EVALUATION OF MODULATORS 

All modulators presented in section 1.2 have been tested with the 
force/torque demand profiles shown in section 1.5. Those profiles have been 
upscaled by different scaling factors until those demands could no longer by 
5 fulfilled because of the limited thrust available. The benefit of using the null 
space search algorithm in the invention (see section 1.4.1) on top of each 
modulator (except the linear optimisation modulator which always yields the 
best possible result) has been analysed. Then, all results are compared with 
respect to the computation time required In order to judge the various 
10 modulators in an overall view. 

Since there Is a lot of data to check, the results are presented in a 
compact form: Because the mass flow scaling (section 1.3) guarantees the 
mass flow constraint to be satisfied and the inventive algorithm presented in 
section 1.4.2 guarantees the satisfaction of the minimum mass flow constraint 

15 on the cost of scaling down the demanded torques and forces, the only 
performance parameter is the demand down scaling factor A<j (1.4.14). From all 
demands of one profile, for every modulator the smallest down scaling factor is 
taken and displayed as performance parameter "100*(1-A<i)%". Clearly, the goal 
is to achieve a loss of 0%. All modulators are compared with respect to loss of 

20 performance and the average computation time required including null space 
search. 

1.6.1 RESULTS 

1.6.1.1 PROFILE #1 

Referring next to Figure 6, there is shown a typical e^ror signal force 
25 demand-thrust characteristic provided for profile #1 . By feeding profile #1 
through the linear pseudo-inverse modulator, the error signal between thrust 
and demand is zero proving that the conrect thrust is generated; this means that 
the demand down scaling factor is one. 

Referring next to Figure 7 there is shown a typical thruster time history as 
30 applied to a linear pseudo-inverse modulator for profile #1 . Note that in the 
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thruster time history, as shown in the Figure, the minimum mass flow constraint 
is satisfied for every thruster because each thruster produces a thrust above the 
"lower bound". 

5 Typical results for other modulators with and without use of the inventive 

null space search are shown in Figure 8. in the results, as shown in Figure 8, 
two scaled up versions of profile #1 have been tested, one with a factor of 2.0 
(solid line) and the other with a factor of 2.5 (dashed line). 

Consider the two diagrams of Figure 8 on the left hand side first (the 
10 lower diagrams are zoom-ins of the upper diagram), showing the result without 
using null space search: Scaling up profile #1 with 2.0, the only modulator 
generating correctly the thrust is the linear optimisation modulator (rectangle- 
symbol), which is also by far the slowest algorithm. The so called generic-type 
extended method (diamond-symbol) comes close (3% loss in demand, this is 
15 the highest loss for all demands in profile #1 with scale factor 2.01), but on the 
other hand is 30 times faster (It needs only 0.014s computation time) The other 
kinds of modulator shown lose 20 % (linear pseudo-inverse (pentagram-symbol) 
and so called generic type simple method (star-symbol)) and 38 % (non-linear 
pseudo inverse method, circle symbol). - Scaling up profile #1 with 2.5 shows 
20 that the linear optimisation method needs to relax the force and torque demand 
as well, because there is not sufficient thrust available to generate this demand. 

Using the inventive null space search to improve the performance (see 
diagrams of Figure 8 on the right hand side) advantageously improves the 
performance of the linear pseudo-inverse and so called generic type simple 
25 method so well that the demand (twice the demand extracted from simulation, 
thus leaving a good safety margin) can be generated by both methods, without 
adding significant computation time. This means that the Inventive null space 
search algorithm conveniently converges very fast. 

The results are encouraging; therefore, the magnetic disturbance torque 
30 has been added (profile #2) to see if the modulators go along with that (higher) 
demand as well. 
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1.6.1.2 PROFILE #2 

Referring next to Figure 9, there are shown some further typical results 
for the various modulators, with and without use of the Inventive null space 
search, thereby generating thrust according to profile #2. 

5 Again, as shown in the Figure, using the inventive null space search 

significantly improves the result of the linear and the non-linear pseudo-inverse 
as well as the so called generic type extended method such that the required 
demand can be generated up to a level of 1.4 times the profile #2. The physical 
limitation is below a scaling of 1.8 times profile #2. 

10 The fastest combination is the non-linear pseudo-inverse method plus 

use of ttie Inventive null space search. 

1.6.1.3 PROFILE #3 

Referring next to Figure 10, there are shown some further typical results 
for the various modulators, with and without use of the inventive null space 
1 5 search, thereby generating thrust according to profile #3. 

This "extreme" profile (it is not generated by STEP-simulation) is a hard 
test for all of the modulators used. 

With a loss of perfonnance of 30%, as shown, the linear pseudo-inverse 
method plus use of the Inventive null space search yields the best result. 

20 1.6.1.4 PROFILE #4 

Refening finally to Figure 11, there are shown some further typical results for 
the modulators, with and without use of the inventive null space search, thereby 
generating thrust according to profile #4. 

Note that the so called generic-type extended method has been 

25 designed particulariy for the case of nonnegative components of the 
force/torque demand and gives by far the best result, with and without use of 
the inventive null space search. Wori<ing out this modulator for all sign 
combinations on the force/torque demand is promising but has not further been 
pursued in this phase of the study. 
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1.7 CONCLUSION ON USE OF MODULATORS IN THE INVENTION 

1) No modulator/null space search combination (apart from the linear 
optimisation modulator which merely serves as a performance reference) 
picks up the total available thrust. 

5 2) To add the inventive null space search after any of the modulators 
significantly improves the result without too much computational effort for 
all modulators (all profiles). This is a definite benefit associated with the 
present invention. 

3) Based on the results, the best and fastest methods are the non-linear 
10 pseudo-inverse method (best for profile #2) and the linear pseudo* 

inverse method (best choice for profile #1) in combination with the 
inventive null space search. They come close to the maximal possible 
performance. For all simulations within STEP, as modulator the linear 
pseudo-inverse method in combination with inventive null space search 
15 has been chosen. 

4) Extreme force/torque demands as in profile #3 can only be achieved with 
a high percentage loss in performance. 

However, it is to be appreciated that all conclusions can hardly be 
generalised and can only be drawn with respect to a certain profile. For a 
20 different profile, the conclusions might be different. Therefore, as the results on 
profile # 4 suggest, the design of the so called generic-type extended modulator 
should be completed, because the analytical formulation allows to derive 
allowable bounds on the components on any force/torque demand for which this 
modulator will work. 

25 Having thus described the present invention by reference to a preferred 

embodiment, it is to be appreciated that the embodiment is in all respects 
exemplary and that modifications and variations are possible without departure 
from the spirit and scope of the invention. For example, the precision of the 
embodiment could possible be improved, if desired, by variation of the number 

30 of thmsters, the only critical requirement being that there are at least two 
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thrusters for use in the Invention. Furtlier, it is to be appreciated tliat in the 
invention the predetenmined constraints couid be simply modified for application 
to other kinds of thruster, ion thrusters for example. 

Furthermore, it Is to be noted that the inventive null space search 
5 algorithm of the Invention is compatible for application to various kinds of 
modulator and can thus be used to improve upon different thruster signals 
proposed by different modulators. 

It is also to be appreciated that tiie thrust demand balancing step could 
be effected in dependence upon tiie modulation step, or alternatively, the 
10 balancing step could be effected wholly independentiy of the modulation step. 

It is to be understood that the number of components on any force/torque 
demand, for example, for conversion Into thrust demands could readily be 
varied to realise the technical effect of the Invention. 

The present invention finds utility for various applications where thrust 
15 demands are required to be balanced, for example balancing thnjst demands in 
spacecraft systems. 



20 



25 
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CLAIMS 



5 
10 

15 2. 

20 3. 
4. 

25 



A method of balancing thrust demands comprising the steps of: 

(a) controllably generating a number of thmst command signals for 
application to a plurality of thrusters; 

(b) comparing the levels of the generated thrust command signals in 
relation to a number of predetemnined constraints; 

(c) identifying which of the thrust demands associated with said 
plurality of thrusters can operate In accordance with said 
predetermined constraints and which of the thmst demands 
associated with said plurality of thrusters can not operate in 
accordance with said predetemnined constraints; and 

(d) balancing the various thrust demands such as to permit each of 
the thmsters to operate in accordance with said predetermined 
constraints. 

A method as claimed in Claim 1, wherein said balancing step comprises 
shifting a predetemiined amount of thrust demand associated with 
particular thrusters operating in accordance with said predetermined 
constraints to one or more of the thrusters which do not operate in 
accordance with said predetermined constraints. 

A method as claimed in Claim 1 or Claim 2. wherein said balancing step 
is effected by application of a predetermined combination of null space 
vectors. 

A method as claimed in Claim 1, 2 or 3, wherein said predetermined 
constraints are selected so that (i) the thrust demand associated with 
each of said plurality of thrusters is more than a predetermined value and 
(ii) the total mass flow associated with said plurality of thrusters is 
constant. 
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5. A method as claimed in any of tlie preceding claims, further comprising 
the step of modulating a number of forces/torques to provide said 
number of thrust command signals. 

6. A method as claimed in Claim 5, wherein said balancing step and said 
5 modulation step are effected independently. 

7. A method as claimed In Claim 5 or 6, wherein said modulation step is 
effected by means of a non-linear pseudo-inverse modulator. 

8. A method substantially as herein described with reference to the 
accompanying drawings. 

10 9. A spacecraft system adapted and anranged to canry out a method as 
claimed in any of the preceding claims. 

10. A system for balancing thrust demands comprising: 

means for controllably generating a number of thrust command signals 
for application to a plurality of thrusters; 
15 means for comparing the levels of the generated thrust command signals 

in relation to a number of predetermined constraints; 

means for identifying which of the thrust demands associated with said 
plurality of thnjsters can operate in accordance with said predetenmined 
constraints and which of the thrust demands associated with said 
20 plurality of thrusters can not operate in accordance with said 

predetermined constraints; and 

means for balancing the various thmst demands such as to permit each 
of the thrusters to operate in accordance with said predetermined 
constraints. 

25 11. A computer program which when loaded into a computer will enable it to 
operate In a system as claimed in claim 10. 
12. A system substantially as herein described with reference to the 
accompanying drawings. 
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Fig.6. 

Error signal force demand-thrust provided for profile #1 
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Thrusters for profile #1 
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